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ABSTRACT 



It is sometimes necessary to know the exact isotopic 
composition of radioactive preparations. Non-radioactive 
isotopes in small amounts may be present in these prepar- 
ations. In this investigation, activation analysis tech- 
niques were applied to carrier free 1-131 to detect and to 
measure the 1-127 present. Although activation analysis 
has been applied extensively to detecting impurities in 
inert materials, it has not been employed generally for 
the detection of isotopes or impurities in radioactive 
preparations. An 1-131 sample was irradiated with thermal 
neutrons at a flux of 5 x 10^neutrons/cm^/sec . The gamma 
spectrum of the irradiated material was determined, and 

the quantity of 1-128 produced was computed. The amount 

18 

of 1-128 found indicated the presence of 5.7 x 10 atoms 

of 1-127 per millicurie of 1-131, or a concentration of 

5 

inactive iodine 1.5 x 10 times that of the radio-iodine. 

The writers are indebted to Professor William W. Hawes 
of the United States Naval Postgraduate School for his pro- 
fessional assistance and his patient encouragement in this 
investigation. 



ii 



TABLE OF CONTENTS 



Section 


Title 


Page 


1 . 


Introduction 


1 


2. 


Neutron Activation 


5 


3. 


Scintillation Spectrometry 


8 


4. 


Analysis 


16 


5. 


Equipment Description 


21 


6. 


Conclusions 


26 


7. 


A Selected Bibliography 


28 






iii 



LIST OF ILLUSTRATIONS 



Figure 

1. 


Decay Schemes of 1-131 and 1-128 


Page 

4 


2. 


Spectral Shape of 1-131 


9 


3. 


Typical Calibration Plot 


11 


4. 


Spectral Shape of 1-128 


13 


5. 


Composite Shape from Irradiated 1-131 


15 


6. 


Sample Calculation 


20 


7. 


Sample Holder 


22 


8. 


Reactor and Console 


23 


9. 


Spectrometer Room 


25 



iv 



1 . INTRODUCTION 






The technique of employing neutron activation for the 
detection of trace element impurities in another material 
began in 1936 with the work of Hevesy and Levi ^lj in which 
a small amount of dysprosium was detected in a sample of 
yttrium. The utility of this method, including its limits 
of detection, was summarized in 1953 by Leddicotte and 
Reynolds J^j . The major advantages of the method as 
shown by these authors are its extreme sensitivity, its 
tendency to be specific because the characteristic radi- 
ations are known for many radio-isotopes, the requirement 
of only small samples, the fact that testing is frequently 
non-destructive in nature, and the accuracy attainable may 
greatly surpass that of other methods. 

Leddicotte and Reynolds report no applications of 
activation analysis techniques for the detection of non- 
radioactive trace elements in a radioactive material. If 
one considers the situation in which the inert trace sub- 
stance is an isotope of the same element as the bulk of the 
radioactive material, it is apparent that normal chemical 
analysis techniques are inapplicable. Activation analysis, 
however, offers some promise for this problem. In this in- 
vestigation, the writers, have sought to determine the feasi- 
bility of applying activation analysis to the detection of 
inert isotopic impurities in radioactive preparations. 
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It may be expected that the relatively intense emiss- 
ions of the bulk of the material will mask those of the acti- 
vated impurity. To a large extent this difficulty can be 
overcome by employing gamma-ray spectrometry. In this way 
the characteristic spectra may facilitate detection. This 
will depend on the spectrometer resolution which becomes im- 
portant when the major gamma emissions of the bulk material 
and the activated isotope are of nearly the same energy. 

1-131 is obtained by chemical separation from fission 
products, and is thus "contaminated" by any other iodine 
isotopes which are produced in fission Ol . When the 
separation is accomplished without adding iodine, the prep- 
aration is referred to as being carrier-free. Isotopes pres- 
ent after moderate cooling tine are inactive 1-127 and the 
long half-life isotope, 1-129, in addition to 1-131 and 1-133 
By radiochemical assay, the purity of carrier-free 1-131 is 
reported to be better than 99 percent 1-131 and less than 
1 percent 1-133. 

If 1-127 is to be detected by activation, it must first 
have sufficiently large thermal neutron cross section for ade 
quate amounts of 1-128 to be formed; and secondly, the major 
gamma emission energy of 1-128 must be sufficiently far from 
the major emissions of 1-131 and 1-132 for the spectrometer 
to resolve it. The cross section of 1-127 for thermal neu- 
trons is about 6 barns (6xl0"^ 4 cm^) which, although small, 
appeared sufficient for activation to be successful with the 
moderate flux available. 
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The major gamma ray emissions of 1-131, 1-132, and 1-128 
are a 0.364 Mev gamma ray, a C.670 Mev gamma ray, and a 
C.455 Mev gamma ray respectively. That is, these particular 
gamma rays are produced after each isotope has decayed by a 
beta process; hence, they represent the decay of excited 
xenon isotopes. (See Figure l) Since these energies differ 
by about 20 percent, and since resolution to within 10 per- 
cent or better can be attained, resolution of the 1-128 
emission as a distinct spectral quantity is assured. 

The above considerations led the authors to select 
carrier-free 1-131 as the radioactive preparation to be 
investigated, and 1-127 as the inert impurity to be detected. 
In brief, the procedure employed consisted of the following 
steps: First, a sample of 1-131 was irradiated in the AGN 

201 Reactor. The sample was then inserted in a Tracerlab 
Stepwise Scanning Spectrometer. The spectrum of the sample 
thus obtained was analized for the presence of the charact- 
eristic 0.455 Mev peak due to 1-128. The observed peak, 
when corrected for decay and for the 1-131 background, 
permitted the amount of 1-127 originally present to be 
calculated . 
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2. NEUTRON ACTIVATION 

The amount of 1-128 produced in an irradiation is given 
by the relationship: 



N 



128 



N 127 °"l27 $ 



X 



( 1 - e 



-Xl28t 



128 



where N 



128 



is the number of 1-128 present after an irradia- 



tion period, t, at a flux level <f> . The terms and X^q 

are the thermal neutron cross section of 1-127 and the decay 
constant of 1-128, respectively. Since Xj^g= 2.77x10“^ rnin*^, 
it can be seen that five half-lives, or 125 minutes, of irrad- 
iation will make the value of this term 97 percent of its 
asymptotic maximum. If the half-life were very long, impract- 
ically long irradiations would be required. Very short half- 
life isotopes on the other hand decay too rapidly for effect- 
ive detection. Leddicotte and Reynolds recommended isotopes 
with half-lives ranging from 10 minutes to several days, as 
being suited for activation analysis. 1-128 with a half-life 
of 25 minutes is within this range. 

The 1-131 preparation was obtained from Oak Ridge Nation- 
al Laboratory, Radio-isotopes Sales Department. The first 
sample of 1-131 consisted of 0.3ml in a glass bottle, and was 
assayed at 34.94mc/ml. It was diluted to a volume of 10ml, 

( lmc/ml ). From this solution 1.8ml were transferred to the 
sample holder for irradiation. The sample holder is illustra- 
ted in Figure 7 and described in Section 5. 
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Each sample was inserted in the reactor via the "glory 
hole" and positioned at the center of the core, where the 
thermal neutron flux is a maximum. All sample irradiations 
were conducted at a power level of 100 milliwatts. The cen- 
tral thermal neutron flux corresponding to this power level 
is rated as 4.5x10^ neutrons/cm^/sec . Following the 125 min- 
ute irradiation, the sample was removed from the reactor and 
transferred to the scintillation spectrometer for analysis. 

The spectroscopic examination of the first sample after 
irradiation resulted in a conclusion that the 1-131 back- 
ground was masking the 1-128, Another sample was prepared by 
concentrating the contents of six vials, each of which con- 
tained approximately 100 microcuries of 1-131 that had been 
decaying for about two years. It was felt that, by lowering 
the background due to 1-131, the presence or absence of 1-127 
could be more clearly established. The sample was irradiated 
and the presence of 1-127 was clearly substantiated. Thus, it 
became necessary to establish a threshold of detection. This 
was done by the activation of samples whose 1-131 activity was 
made progressively less either by dilution or by decay. 

Other irradiations were performed to generate data neces- 
sary for the quantitative analysis. An empty sample holder 
was irradiated and analized; no appreciable activity was ob- 
served. A solution containing a known amount of 1-127 was ir- 
radiated in order to generate the 1-128 spectral shape. This 
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sample was also used to obtain a calibration between the ob- 
served 1-128 activity in the 1-131 samples and the amount of 
1-127 present. 



•n 
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3. 



SCINTILLATION SPECTROMETRY. 



A recent report on scintillation spectrometry by R.L. 
Heath [" 3 J proved to be most applicable for this investigat- 
ion. In this, a catalogue of gamma ray spectra, the concept 
of spectral shapes is clearly explained, and the reader is 
referred to this article for details. A given radioactive 
nuclide when examined with a spectrometer, will yield an 
intensity distribution as a function of energy for the 
emitted gammas, their Compton smears, and backscattered rays. 
If the backscattering contribution is reduced greatly by the 
shielding arrangement, or eliminated mathematically, then the 
intensity distribution curve as a function of energy is call- 
ed the characteristic "spectral shape" of the nuclide. (See 
Figure 2). If more than one nuclide is present in a sample, 
the obtained envelope is a composite of the appropriately 
weighted individual spectral shapes. 

Intensity-versus-energy measurements were made for sev- 
eral types of samples. The sample types were: background, an 
irradiated empty sample holder, a sample of irriadiated 1-127, 
a sample of non-irradiated 1-131, and samples of irradiated 
1-131 with various specific activities and decay times. 

■r 

A number of spectrometer gain, voltage, discriminator, 
and bias settings were examined. The most linear and most 
stable response with suitable resolution was obtained when 
these settings were such as to give a range of 700 Kev for 
the 50 channels, or a 14 Kev increment per channel. A Sola 



8 



RELATIVE INTENSITY 




FIGURE 2 
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transformer was interposed between the soectrometer and the 
line to insure good voltage stability. All chassis compon- 
ents were kept on continuously to insure maximum thermal sta- 
bility . 

The resolution time of the spectrometer is such that 
coincidence losses become significant for total count rate in 
excess of 4000 counts per second. During the examination of 
the higher specific activity, irradiated, 1-131 samples, this 
limiting value was closely approached. 

The response of the amplifier is rated as linear within 
1 per cent. Many calibrations at various spectrometer set- 
tings were performed using Pa-234, Zn-65, Co-60, Cs-137, and 
1-131 as calibrating sources. Except for excessively high 
total count rates, as mentioned above, the linearity was ob- 
served to be within the specification. For the final experi- 
mental runs, Cs-137 and 1-131 were used exclusively as the 
calibration sources. (See Figure 3) 

A 50 minute per channel background examination plus 
numerous other shorter background checks showed that the 
average background count rate in each of the first 10 chan- 
nels was about 20 counts per minute. This value gradually 
diminished to essentially zero in the higher channels. Since 
during the major sample runs the channel count rates were in 
the thousands, background effects were not included. 

The irradiated empty sample holder, when counted, showed 
a slight amount of activity, about 2 to 4 counts per minute 
above background. No photopeaks were observed. Hence, the 
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FIGURE 3 



contribution of the sample holder may be ignored. The sam- 
ples of irradiated 1-127 and non-irradia ted 1-131, when ex- 
amined and corrected for decay, gave the spectral shapes 
shown for 1-128 and 1-131 respectively. (See Figures 2 and 4) 
These shapes were used in the analysis of the irradiated 1-131 
samples for the presence of 1-128. Using the spectral shape 
subtraction approach as described by Heath, the effective peak- 
to-total ratio for the 0.455 Mev emission of 1-128 was obtain- 
ed. This peak-to-total ratio for the experimental configur- 
ation employed was 0.185. The subtraction scheme was applied 
to the irradiated 1-131 samples to eliminate the 1-131 con- 
tribution to the composite spectral shape. 

The peak-to-total ratio is the ratio of the area under 
the photopeak to the area under the entire spectral shape, 
provided one has a monoenergetic source. In the event there 
are two or more emission energies in the source, one must use 
a subtraction scheme. This approach involves fitting an idea- 
lized monoenergetic spectral shape, photopeak to photopeak, on 
the highest energy photopeak in the sample spectrum. Once 
fitted, this spectral shape is subtracted from that of the 
sample, leaving the contribution of lesser energy peaks only. 
Successive applications of this approach permit the composite 
spectral shape to be resolved into its components. The peak- 
to-total ratio for each component may then be determined. 

To determine the peak-to-total ratio for a monoenergetic 
emission or for a resolved component, one fits a Gaussian 
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RELATIVE INTENSITY 




CHANNEL 
FIGURE 4 
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distribution to the high energy side of the photopeak. It is 
assumed that the area under the Gaussian distribution is the 
same as that area under the experimental photopeak. The area 
under the balance of the spectral shape may be determined by 
any appropriate method. The relationship of area to maximum 
ordinate enables one to measure, in effect, the area under 
the photopeak by measuring its height. If the peak-to-total 
ratio has previously been determined, then one can immediate- 
ly calculate the number of events of the particular energy 
that occurred in the scintillation crystal. 

The emissions of 1-128 of energies higher than 0.455 Mev 
contribute only 11 percent of the total gamma emission. 
Therefore, by using the spectral shape subtraction technique 
to adjust the 1-128 shape to that of its equivalent monoener- 
getic 0.455 Mev spectral shape, little accuracy is lost. 

Small percentage errors in fitting these higher energy peaks 
prior to subtraction are reduced nearly ten-fold in their 
effect on the equivalent 0.455 Mev spectral shape. In view 
of these considerations it is felt that the value of 0.185 for 

the peak-to-total ratio is accurate to within 1 percent. An 

•v 

example of the irradiated 1-131 preparation with composite spec- 
tral shapes of 1-131 and 1-128 is shown in Figure 5. 
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FIGURE 5 
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4. 



ANALYSIS. 



In the last section the concept of peak-to-total ratio 
was discussed. The area under the photopeak was assumed to 
be that under a Gaussian distribution fitted to the high en- 
ergy side of the photopeak and having the same value for its 
maximum ordinate. Thus, the area under the photopeak is 
2.56 times the value of the maximum ordinate. In terms of 
the emission rate under the photopeak R^, and the total 
emission rate R^., the peak-to-total ratio p is : 



P 



Rp 

R t 



If there were no absorption losses, and if the detection 
efficiency were 100 percent, the absolute emission rate R 
would be equal to the total emission rate R^. The methods 
employed to obtain these corrections will be discussed in 
later paragraphs. The resultant equation for the determin- 
ation of the absolute emission rate R is : 



R4 



Rr 



Rr 



R = 



A e. 



In 

Rt 



Ae H 



p e t A 



where A is the absorption correction, and e^ is the detection 
efficiency. For given emission, p,A, and e^ can be determin- 
ed independently and are independent of the intensity of the 
emission within the normal range of operation. Hence the 
absolute emission rate for any sample can be calculated as 



soon as the maximum ordinate of the photopeak has been deter- 
mined . 



The absorption correction was computed by using the 
gamma-ray absorption coefficient of water for the sample sol- 
ution, and an average of water and carbon gamma-ray absorption 
coefficients for the plastic of the sample holder. The value 
obtained was A = 0.935. 

The total absolute detection efficiency e^ was interpol- 
ated from disk source efficiency calculations given by R.L. 
Heath. . The sample holder volume was partitioned into thin 
layers, or disks, The disk source efficiency for each layer 
was then determined. The mean of the efficiencies thus deter- 
mined was e^ = 0.0256. 

The 0.455 Mev emission of 1-128 accounts for 17.2 per- 



0.455 Mev emission is related to the number of 1-128 atoms 
present, Nj^g, in the following manner : 



The number of 1-128 atoms was always obtained from a decay 
corrected value of R corresponding to the time of termina- 
tion of irradiation. N^ 2 g is related to the original numbe] 
of 1-127 atoms present, N i27 , b Y the relationship, as shown 
on page 5. 



cent of the total decay. The decay constant \^ 2 g of 

1-128 is 2.77 x 10“^ min"!. The absolute emission rate R for 




R 



0.172 \ 



N 128 = 209 • 9 R 



128 
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The examination of the irradiated 1-131 for the pre- 
sence of 1-128 was accomplished by taking measurements only 

\ 

in the channels in the immediate vicinity of the 0.455 Mev 
photopeak, and recycling over these same channels as rapidly 
as possible during the decay. All such measured count rates 
were decay corrected to the time of the end of irradiation. 
After this correction all count rates for each respective- 
channel were averaged. These averaged values were then plot- 
ted and the spectral shape subtraction technique applied to 
obtain the 1-128 contribution to the composite spectrum. 

For the 1-131. sample taken from the two year old prep- 
aration, no accurate value of the amount of 1-131 originally 
present was known. Therefore, any calculation of. the amount 
of 1-127 in this sample is not relatable to the amount of 
1-127 per millicurie of 1-131. 

Another sample of 1-131 had been decaying for 36 days. 

It consisted of 5 x 10“^ ml of stock solution diluted to fill 
the sample holder. The stock solution was assayed at Oak 
Ridge as being 29 mc/ml of 1-131. Hence at the time of ex- 
amination the specific activity of the stock solution was 
1.35 mc/ml. The spectral shape obtained indicated the pre- 
sence of 1-128, but was not sufficiently well defined for re- 
liable evaluation of 1-127 content. This, then, constitutes 
the detection limiting activity. 

A third sample was prepared from an 1-131 solution that 
had decayed for 81 days. A local assay of this solution agreed 

precisely with the Oak Ridge assay of 35 mc/ml for 0.3 ml. 
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The 0.3 ml was diluted to 10 ml, of which 6 ml were taken, 
concentrated, and put into a sample holder. This sample was 
irradiated, and then examined as described above. The decay 
corrected photopeak was 42 ± 2 cpm/unit energy interval, so 
that Rp = 108 ± 5 cpm. This represents 8.6 x 10^ atoms of 
1-127 in the sample, or 18.1 ± 0.9 mg. A sample of 14 mg of 
1-127 was irradiated under the same circumstances as the 
1-131 sample. The decay corrected photopeak height was 81.5 
cpm/unit energy interval. Using this direct calibration the 
amount of 1-127 in the 1-131 sample was 7.49 mg as compared 
to the calculated value of 18,1 mg. The calculated value is 
2.4 times the direct comparison value. It is probable that 
the discrepancy is due to error in the assumed value of the 
flux. See Figure 6 for a sample calculation. 






SAMPLE CALCULATION 



VALUES OF CONSTANTS: 

p = 0.185 e t = 0.0256 A = 0.935 

= 4.5 x lC^n/cm 2 /sec ^128 = 2.77' x I0~ 3 min”J- 

6 ^ 27 = 6.3 x 10~ 24 cm 3 = 4.62 x 10" 4 sec“l 

OBSERVED HEIGHT OF 1-128 PHOTOPEAK : 

42 cpm/unit energy interval. 

AREA UNDER 1-128 PHOTOPEAK : 

Rp = 2.56 x 42 = 108 cpm 

ABSOLUTE EMISSION RATE FOR 0.455 MEV GAMMA OF 1-128 : 

Rp A 

R = - = 2.42 x 10 4 cpm 

p e t A 

NUMBER OF 1-128 ATOMS AT END OF IRRADIATION : 



NUMBER 



NUMBER 



N, OQ = = 209.9 R = 5.09 x 

128 °- 172 X 128 
OF 1-127 ATOMS IN SAMPLE : 

= X 128 N 128 x 1 _ 

127 127 ' ( 1 - e” X l28 t 

8.55 x 10 ^ atoms 

OF MG OF 1-127 IN SAMPLE : 

N 127 X A 127 ^ 8.55 x 10 19 X 127 x 10 3 

6.02 x 10 23 6.02 x 10 23 

FIGURE 6 



lp^ atoms 

A 



= 18.1 mg 
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5. EQUIPMENT DESCRIPTION. 

The Aerojet-General Nucleonics Model 201 Research 
Reactor has been described ["bj . The reactor, when operating 
at a rated power level of 100 milliwatts, has a central core 
thermal neutron flux 4.5x10^ neutrons/cm^/sec , a flux level 
comparable to that attainable with moderate intensity Ra-Be 
or Po-Be sources. 

The "glory hole"of the AGN 201 reactor is an aluminum 
tube, approximately one inch in diameter, which passes 
horizontally through the center of the reactor. Two plug 
rods are normally inserted in the "glory hole". These rods 
are so constructed that the materials of their inner ends 
match the reactor core, reflector, and gamma ray shield for 
that volume not occupied by samples being irradiated. The 
general arrangement of the reactor and its control console 
are shown in Figure 8. 

The sample holders employed were designed to pass 
freely through the "glory hole". Each holder was assembled 
and fused with glacial acetic acid. The dimensions and 
design of the sample holders are shown in Figure 7. As a 
further safeguard against leakage the sample holder was 
sprayed with an acrylon plastic sealer. 

A Tracerlab Stepwise Scanning Spectrometer, RLP-6 III, 
was employed in this investigation. Its major components 
are: a Nal (Tl) scintillation detector, a power supply, a 
high voltage supply, a pulse amplifier and shaper, a thres- 
hold discriminator, a stepwise-scanning single window 

21 
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This photograph shows the AGN 201 Reactor and 
its control console. Lt. Gminder is inserting the 
plug rod into the "glory hole". Major Tyler is at 
the control console. Professor Hawes is supervising 
the operation. 



REACTOR AND CONSOLE 



FIGURE 8 
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discriminator, and an automatic read-out scaler and printer. 
The detector head consists of a lJ^"xl" Nal(Tl) scintillation 
crystal coupled to Dumont Type 6292 photomultiplier tube, 
both components being housed in a light-proof metal casing. 
This housing is coupled to a sample stage, and the whole 
assembly is encased by a three inch thick cylindrical lead 
shield. The high voltage supply provides a well regulated 
voltage from 500 to 5000 volts in 10 volt increments. Its 
stability is rated as less than 0.02 percent drift per day, 
and its regulation is such that a. 1 percent change in line 
voltage will cause no more than 0.2 volts change in the 
output high voltage. The pulse amplifier is linear to with- 
in 1 percent and has a rise time of 0.25 p.sec. The stepwise 
scanning element moves a fixed window, 5 volts in width, 
through 50 steps, or channels, covering a range of 0 to 
255 volts. The automatic read-out scaler and printer work 
in conjunction with the scanning element so that as inform- 
ation obtained for a given channel is recorded, the scanning 
element advances to the next higher channel. Upon completion 
of each recording cycle, counting in the next higher channel 
commences. The read-out information consists of the channel 
number, the number of counts occurring in this channel during 
a selected interval of time, and the time interval in minutes. 
The spectrometer arrangement is shown in Figure 9. 
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This photograph shows the spectrometer 
From left to right, the spectrometer system 
the re. id-out typewriter, the RLP-6 cabinet, 
shield wit.', the detector unit inside. 



arrangement . 
consists of 
and the lead 
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CONCLUSIONS. 



This investigation has shown that an inert isotope, 
1-127, is present in carrier-free 1-131, For the specific 
sample that was quantitatively measured this 1-127 content 
was found to be 1.2 mg of 1-127 per me of 1-131. This inves- 
tigation has demonstrated that studies of this nature can be 
successfully pursued with low intensity neutron sources pro- 
vided that gamma ray spectrometry is employed for radiation 
analysis. This conclusion is supported by recent studies at 
the U.S. Naval Hospital at Bethesda, Maryland . jioj . Workers 
there have sought to determine the carrier-free status of 
radioisotopes by neutron activation. Their investigation 
was specifically applied to the detection of Cu-63 in Cu-64, 
detection being indicated by the increase in activity. Such 
increases were reported, indicating that Cu-63 is present in 
carrier-free Cu-64. The techniques utilized here, however, 
should greatly improve the sensitivity of that determination. 
Since the sensitivity of this analysis is directly proport- 
ional to the neutron flux level, one may conclude that, by 
increasing the flux level from 10^ to 10 ^ neutrons/cm^/sec , 
a sensitivity of 0.018 pg may be achieved. This is comparabl 
to the 0.002 \ig sensitivity quoted by Leddicotte and Reynolds 
with inactive materials for 1-128, using a flux of this mag- 



nitude 



.[ 2 . 12 ], 



The successful application of neutron activation 

analysis to a radioactive preparation is dependent mainly on 

two factors. First of all, the combination of the thermal 
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neutron cross section of the nuclide sought and the avail- 
able neutron flux level must be such as to assure that a de- 
tectable activity will be produced. Secondly, the emission 
of the radioactive preparation, its decay products, and the 
activated products other than that sought must not be suffi- 
ciently close, energy-wise, to that of the sought activated 
nuclide so as to mask its presence. This second factor is 
cited by Wahl and Bonner as being one of the two major limit- 
ations of activation analysis. M. 
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